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Abstract 

We are concerned with the discretization of a solution of a Forward-Backward stochastic 
differential equation (FBSDE) with a jump process depending on the Brownian motion. In 
this paper, we study the cases of Lipschitz generators and the generators with a quadratic 
growth w.r.t. the variable a. We propose a recursive scheme based on a general existence 
result given in the companion paper [15] and we study the error induced by the time 
discretization. We prove the convergence of the scheme when the number of time steps n 
goes to infinity. Our approach allows to get a convergence rate similar to that of schemes 
of Brownian FBSDEs. 
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1 Introduction 


In this paper, we study a discrete-time approximation for the solution of a forward-backward 
stochastic differential equation (FBSDE) with a jump of the form 


Xt = X + f b{s,Xs)ds+ f a{s,Xs)dWs+ f P{s,Xs-)dHs , 

Jo Jo Jo 

Yt = g{XT) + j^ fis,X,,Y„Z„Us)ds- Z.dW^- U,dHs , 
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where = !,-<* and r is a jump time, which can represent a default time in credit risk 
or counterparty risk. Such equations naturally appear in finance, see for example Bielecki and 
Jeanblanc [5], Liin and Quenez [T5], Peng and Xu [50], Ankirchner et al. [T] for an application to 
exponential utility maximization problem and Kharroubi and Lim m for the hedging problem 
in a complete market. The approximation of such equation is therefore of important interest 
for practical applications in finance. In this paper, we study the case where the generator / is 
Lipschitz or with quadratic growth w.r.t. Z. 

In the literature, the problem of discretization of FBSDEs with Lipschitz generator has been 
widely studied in the Brownian framework, i.e. no jump, see e.g. [I1I71I3I1I11IS]. More 
recently, the case of quadratic generators w.r.t. Z has been considered by Imkeller et al. [8] and 
Richou m- For Lipschitz generators, the discrete-time approximation of FBSDEs with jumps is 
studied by Bouchard and Elie |1] in the case of Poissonian jumps independent of the Brownian 
motion. Their approach is based on a regularity result for the process Z, which is given by 
Malliavin calculus tools. This regularity result for the process Z was first proved by Zhang [55] 
in a Brownian framework to provide a convergence rate for the discrete-time approximation of 
EBSDEs. The use of Malliavin calculus to prove regularity on Z is possible in [4] since the 
authors suppose that the Brownian motion is independent of the jump measure. 

In our case, we only assume that the random jump time r admits a conditional density given 
W, which is assumed to be absolutely continuous w.r.t. the Lebesgue measure. In particular, 
we do not specify a particular law for r and we do not assume that r is independent of W as 
for the case of a Poisson random measure. 

To the best of our knowledge, no Malliavin calculus theory has been set for such a framework. 
Thus, the method used in [1] fails to provide a convergence rate for the approximation in this 
context. 

We therefore follow another approach, which consists in using the decomposition result given 
in the companion paper m to write the solution of a PBSDE with a jump as a combination 
of solutions to a recursive system of FBSDEs without jump. We then prove a regularity result 
on the Z components of Brownian BSDEs coming from the decomposition of the BSDE with a 
jump. This regularity result allows to get a rate for the convergence of the discrete-time schemes 
for these BSDEs as in [55] or [4] for the Lipschitz case and [21] for the quadratic case. 

Finally, we recombine the approximations of the solutions to recursive system of Brownian 
FBSDEs to get a discretization of the solution to the FBSDE with a jump. 

We notice that our approach also allows to weaken the assumption on the forward jump 
coefficient in the Lipschitz case. More precisely, we only assume that /3 is Lipschitz continuous, 
unlike [4] supposing that /3 is regular and the matrix -I- V/3 is elliptic. 

As said above, this kind of FBSDEs with a jump appears in finance. The general assumptions 
made on the jump time r allow to modelize general phenomena as a firm default or simpler 
as a jump of an asset that can be seen as contagion from the default of another firm on the 
market, see e.g. m for some examples. In particular, the approximation of these FBSDEs has 
its own interest, since it provides approximations of optimal gains and strategies of the studied 
investment problems. 

We choose to present our results in the case of a single jump and a one-dimensional Brownian 
motion for the sake of simplicity. We notice that they can easily be extended to the case of a 
d-dimensional Brownian motion and multiple jumps with eventually random marks, as in |15j . 
taking values in a finite space. 

The paper is organized as follows. The next section presents the framework of progressive 
enlargement of a Brownian filtration by a random jump, and the well posedness of FBSDEs in 
this context. In Section 3, we present the discrete-time schemes for the forward and backward 
solutions based on the decomposition given in the previous section. Finally, in Section 4, we 
study the convergence rate of the scheme for the forward solution. In Sections 5 and 6, we study 
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the convergence rate of the scheme for the backward solution for the Lipschitz case respectively 
for the quadratic case. 

2 Preliminaries 

2.1 Notation 

Throughout this paper, we let (n,C/,P) a complete probability space on which is defined a 
standard one dimensional Brownian motion W. We denote F := the natural filtration 

of W augmented by all the P-null sets. We also consider on this space a random time r, i.e. 
a nonnegative J^-measurable random variable, and we denote classically the associated jump 
process by H which is given by 


Ht := lr<i , t > 0 . 

We denote by D := {'Dt)t>o the smallest right-continuous filtration for which r is a stopping 
time. The global information is then defined by the progressive enlargement G := (0t)t>o 
the initial filtration where 

Gt ■= Pi V 'Dt+e^ 

e>0 

for all t > 0. This kind of enlargement was introduced by Jacod, Jeulin and Yor in the 80s (see 
e.g. m, [H] and M)- We introduce some notations used throughout the paper 

- 7^(F) (resp. 'P(G)) is the cr-algebra of F (resp. G)-predictable measurable subsets of 
G X R+, i.e. the cr-algebra generated by the left-continuous F (resp. G)-adapted processes, 

- 'PAf(F) (resp. 'PAf(G)) is the tr-algebra of F (resp. G)-progressively measurable subsets 
of G X 1R+. 

We shall make, throughout the sequel, the standing assumption in the progressive enlargement 
of filtrations known as density assumption (see e.g. [laiiiiiii]). 

(DH) There exists a positive and bounded 7^(F) (g) S(M+(-measurable process 7 such that 

F[t G d0 I J-t] = ')t{0)dO , t > 0 . 


Using Proposition 2.1 in m we get that 
intensity. 


(DH) ensures that the process H admits an 


Proposition 2.1. The process H admits a compensator of the form Xtdt, where the process A 
is defined by 


Xt 


Itjt) 

P[t > t \ Tt] 


1 


t<T 5 


t>0 . 


We impose the following assumption to the process A. 

(HBI) The process A is bounded. 

We also introduce the martingale invariance assumption known as the (H)-hypothesis. 

(H) Any F-martingale remains a G-martingale. 

We now introduce the following spaces, where a,b G IR_|_ with a < b, and T < 00 is the terminal 
time. 
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- iS|^[a,6] (resp. 5“[a,6]) is the set of VM{G) (resp. 7^A1(F))-ineasurable processes 
(Xt)te[a,b] essentially bounded 

II^IU“[a.b] := esssuplYtl < cx) . 


- iSg[a, 6] (resp. <Sp[a, &]), with p > 2, is the set of VM{G) (resp. "PAl(F))-measurable 
processes (Ft)tG[a.b] such that 


ll^lUpfa.b] 


1 



- H^[a, b] (resp. H^[a, 6]), with p > 2, is the set of V{G) (resp. 7’(F))-measurable processes 
{Zt)te[a,b] such that 


\HP{a,b] 



< OO . 


- i^(A) is the set of 'P(G)-measurable processes {Ut)t^[o,T] such that 


WuWlh,) 




< OO . 


2.2 Forward-Backward SDE with a jump 

Given measurable functions b : [0, T] x R —)■ K., cr : [0, T] x R ^ R, /3 : [0, T] x R —)> R, p : R —R 
and / : [0, T] x R x R x R x R —>■ R, and an initial condition a; G R, we study the discrete¬ 
time approximation of the solution {X, Y, Z, U) in <Sg[0, T] x *S^[0, T] x T] x L^(A) to the 

following forward-backward stochastic differential equation 


Yt 


X-\- fb{s,Xs)ds+ f a{s, Xs)dWs + fl3{s,Xs-)dHs , 0<t<T. 

Jo Jo Jo 

g{XT)+ J^f{s,Xs,Ys,Zs,{l-Hs)Us)ds 


- [ ZsdWs - [ UsdHs , 0 < t < T , 


( 2 . 1 ) 


( 2 . 2 ) 


when the generator of the BSDE is Lipschitz or has a quadratic growth w.r.t. Z. 


Remark 2.1. In the BSDE (12.21) . the jump component U of the unknown {Y,Z,U) appears 
in the generator / with the additional multiplicative term 1 — H. This ensures the equation to 
be well posed in 5g^[0,T] x Hq[0,T] x T^(A). Indeed, the component U lives in T^(A), thus 
its value on (r A T, T] is not defined since the intensity A vanishes on (r A T, T]. We therefore 
introduce the term 1 — R to kill the value of U on (r A T, T] and hence to avoid making the 
equation depending on it. 


We first prove that the decoupled system admits a solution. To this end, we 

introduce several assumptions on the coefficients 6 , ct, /3, g and /. We consider the following 
assumption for the forward coefficients. 

(HF) There exists a constant K such that the functions b, a and /? satisfy 


|&(t,0)|-h |cr(t,0)|-h |/3(t,0)| < K, 
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and 


\b{t^x) — b{t^x')\ + \a{t,x) — a{t,x')\ + \l3{t,x) — I3{t,x')\ < K\x — x\^ 

for all (f, X, x') € [0, r] X R X M. 

For the backward coefficients g and /, we impose the following assumptions for the Lipschitz 
case. 

(HBL) There exists a constant K such that the functions g and / satisfy 

|/(t,a:,0,0,0)| + |5(a;)| < K, 

and 

\f{t,x,y,z,u) - f{t,x,y',z',u')\ < K(\y - y'\ + \z - z'\ + \u - u'\) , 

for all (t, X, y, y', z, z', u, u') S [0, T] x R x x x R^. 

For the backward coefficients g and /, we consider the following assumptions for the quadratic 
case. 

(HBQ) 

- There exist three constants Mg, Kg and Kg such that the functions g and / satisfy 

|g(x)| < Mg , 

\g[x)-g[x)\ < Kg\x-x\, 

\f(t,x,y,z,u) - f{t,y',z,u)\ < Kg\y-y'\, 

\f{t,x,y,z,u)\ < Kq{l + \y\ + \z\^ +\u\) , 

for all {t, X, x', y, y', z, u) e [0, T] x R^ x R^ x R x R. 

- For any i? > 0 there exists a function mc;^ such that limE_>o mc^^{e) = 0 and 

\f{t,x,y,z,u-y) - f{t,x,y',z',u-y')\ < 

for all {t,x,y,y',z,z',u) € [0,T] x R x R^ x R^ x R s.t. jyl, \z\, \y'\, \z'\ < R and 
\y-y'\ + \z-z’\<e. 

- f{t,u) = f{t, 0) for all u G R and all t G {t AT, T], 

- The function f(t,x,y, .,u) is convexe (or concave) uniformly in {t,x,y,u) G [0,T] x R x 
R X R. 


In the sequel K denotes a generic constant appearing in (HBL), (HBQ) and (HF) and which 
may vary from line to line. 

In the purpose to prove the existence of a solution to the FBSDE we follow the 

decomposition approach initiated by m and for that we introduce the recursive system of 
FBSDEs associated with (CT)-(I^. 

• Find {X\0),Y\e),Z^{0)) G <S^[0,r] X <S^[6»,r] X H^[0,T] such that 

Xli0) = x + [b{s,Xl{0))ds+ [a{s,Xl{0))dW, + l3{0,Xl.{0))lg<t , 0<t<T, (2.3) 
Jo Jo 

Y,\0) = g{X^{0))+lJ{s,Xl{0),Y,\0),Zl{0),O)ds-lzl{0)dWs, 0<t<T, (2.4) 
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for all e e [0,T]. 


Find {X°,Y°,Z°) G 5|[0,T] x 5^[0,T] x H^[0,T] such that 

X° = x + [ b{s,X°)ds + [ a{s,X°)dWs , 0 < t < T , 

Jo Jo 

= g{Xi^) + r f{s,X°,Y°,ZlY}{s)-Y°)ds- r Z^Ws. 0<t<T. 
Jt Jt 


(2.5) 

( 2 . 6 ) 


Then, the link between the FBSDE (l2.1F(j2.2l) and the recursive system of FBSDEs (I2.3p - (j2.4p 
and (I2.5I1 - (I2.6I) is given by the following result. 


Theorem 2.1. Assume that (DH), (HBI), (H), (HF) and (HBL) or (HBQ) hold true. 
Then, the FBSDE (I2.1II - (I2.2II admits a unique solution {X,Y, Z,U) G 5q[0,T] x iS^[0,r] x 
Hg[0,T] X L'^{X) given by 

' Xt = X?lt<r + Xl{T)lr<t , 

Yt = y,°w + y/(T)i.<t, 

Zt = Z^'i.t<T F Zl{T)lTct , 

^ Ut = {Y,\t)-Y°)lt<r, 

where {X^{9),Y^{0), Z^{6)) is the unique solution to the FBSDE (I2.3I) - (I2.4I) in iSp[0,T] x 
iS]f“[0,T] X F[f[9,T], for 9 G [0,r], and (X^,Y^, Z^) is the unique solution to the FBSDE 
dZSI-JMl) zn>s^[o,r] x5|“[o,r] xi^2[o,^]. 


Proof. 

Step 1. Solution to m under (HF). 

Under (HF) there exist unique processes G iSp[0,T] satisfying (12.51) . and X^{9) G <Sp[0,r] 
satisfying (12.31) for all 9 G [0,r] such that X^ is VM{¥) (g) S(]R+)-measurable. Then, from the 
definition of Ff, we check that the process X defined by 


X* = X^^tt<r + Xl{T)tt>r , 


( 2 . 8 ) 


satisfies m- We now check that X G 5g[0,T]. We first notice that from (HF), there exists a 
constant K such that 


E 


sup 
LtGlO.T] 


< K . 


(2.9) 


Then, from the definition of X^ and X^, we have for all t G [9,T] 

sup \Xl{9)\^ < k(\X^\^ + \l3i9,X°)\^ + [ \b{u,Xl{9))\^du+ sup [ aiu,Xii9))dW^ 


Using (HF) and BDG-inequality, we get 


E 


sup 

se[ 6 (,t] 


|^i(0)| 


< K 





Applying Gronwall’s lemma, we get 


sup ||X1(0)|L,.„„, < K. (2.10) 

ee[o.T] 

Combining (12.81) . (12.91) and (12.101) . we get that X G 5Q[0,r]. Moreover still using (HF) we get 
the uniqueness of a solution to (12.11) in 5,g[0,T]. 
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step 2. Solution to under (DH), (HBI), (H) and (HBL). 

To follow the decomposition approach initiated by the authors in we need the generator 
to be predictable. To this end, we notice that in the BSDE (I2.2L we can replace the generator 
{t,y,z,u) f{t,Xt,y,z,{l-Ht)u) by the predictable map {t,y,z,u) f{t,Xt-,y,z,{l- 

Using the decomposition (j2.8|) . we are able to write explicitly the decompositions of the 
C/T-measurable random variable giXr) and the 'P{G) 0 S(R) 0 S(R) 0 S(]R)-measurable map 
(w, t, y, z, u) fit, Xt- (w), y, z, 14(1 - Ht- (w))) given by Lemma 2.1 in [15] 

giXr) = giX°)lT<r + giX^iT))lT>r , 
f{t,Xt-,y,z,{l-Ht-)u) = fit,y,z,u)lt<r + fit,y,z,u,T)lt>r, 

with/°(t,y,z,ii) = fit,X^,y,z,u) and f^it,y, z,u, 9) = /(t, (0), y, z, 0), for all it,y,z,u,6) € 

[0, T] X K X K X R X R+. 

Suppose now that (DH), (HBI), (H) and (HBL) hold true. Then, from Theorem C.l 
in [15], the BSDE (12.41) admits a P(F) 0 S([0, r])-measurable solution iY^,Z^) and the BSDE 
(|2.6p admits a solution (U*^,0°). Using Proposition 2.1 in we obtain 

r'(0)ll5=»[e.T] + ll^'Wllil=[e.T] < K, 

for all 6 S [ 0 , T], and 

11^ ll5“[0,T] + 11^ IIh2[o,T] — ^ ■ 

We can then apply Theorem 3.1 in [HI and we get the existence of a solution to (ESI) in 
5g°[0,r]xiL2[o,T]xL2(A). 

Let (F, Z, U) and {¥', Z', U') be two solutions to (12.21) in 5^[0, T] x T] x L^(A). Since 
fit, X, y, z, (1 — Ht)u) = fit, X, y, z, 0) for all t 6 (r A T, T] and A vanishes on (r A T, T], we 
can assume w.l.o.g. that Ut = C/( = 0 for t G (r A T,T]. Then, from (DH), (HBI), (H) and 
(HBL), we can apply Theorem 4.1 in m and we get that Y <Y'. Since Y and Y' play the 
same role, we obtain Y = Y'. Identifying the pure jump parts of Y and Y' gives U = U'. 
Finally, identifying the unbounded variation gives Z = Z'. 

Case 3. Solution to (12.21) under (DH), (HBI), (H) and (HBQ). 

The existence of a solution {Y,Z,U) G 5Q[0,r] x Lq[0,T] x L^(A) is a direct consequence of 
Proposition 3.1 in [15]. We then notice that from the definition of H we have /(t, x, y, z,«(! — 
Ht)) = fit,x,y,z,0) for all t G (r Ar,T]. This property and (DH), (HBI), (H) and (HBQ) 
allow to apply Theorem 4.2 in which gives the uniqueness of a solution of ( 1 ^ . □ 

Throughout the sequel, we give an approximation of the solution to the FBSDE (I2.ip - (I2.2I) 
by studying the approximation of the solutions to the recursive system of FBSDEs (I2.3D - (I2.4D 
and (I2.5I) - (I2.6I) . 

3 Discrete-time scheme for the FBSDE 

In this section, we introduce a discrete-time approximation of the solution (X, Y, Z, U) to the 
FBSDE (I2.1D - (I2.2I) based on its decomposition given by Theorem 12.II 

Throughout the sequel, we consider a discretization grid tt := {tg,... ,tn} of [0,T] with 
0 = to < ti <...< tn = T. For t G [0, T], we denote by 7 r(t) the largest element of tt smaller 
than t 

nit) := max { ti , i = 0 ,... ,n \ ti < t } . 
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We also denote by |7r| the mesh of tt 

|7r| := max { tj+i - , z = 0,... ,n - 1 } , 

that we suppose satisfying |7r| < 1, and by (resp. the increment of W (resp. the 

difference) between ti and AW^ := — Wt^_^ (resp. Atf := ti — ti-i), for 1 < z < n. 


3.1 Discrete-time scheme for X 


We introduce an approximation of the process X based on the discretization of the processes 
X° and X\ 

• Euler scheme for We consider the scheme X^’'^ defined by 


\^0,7r 

■^to 


= X 


X 


0,7r _ 


= At- + a(t,_i, X'i::^)XW: , 1 < z < n 


r0,7r 


0,7r 


• Euler scheme for X^. Since the process X^ depends on two parameters t and 9, we introduce 
a discretization of in these two variables. We then consider the following scheme 

( = x + /3{to,x)l^(^e)=o , 

I Xl;-{7r{e)) = Xt^(7r(0))+6(t._i,Xt"(^(0)))At-+a(t._i,At\(^(0)))AW- (3.1) 

[ + P{ti_i,Xl’^^{TT{9)))lt,=^(e) , l<i<n, O<0<T. 


We are now able to provide an approximation of the process X solution to the FSDE (12.11) . We 
consider the scheme X- defined by 


X- = X°(-W+X^(-(^(r))li>. , 0<t<T. (3.2) 

We shall denote by {X°’-}o<i<n (resp. {Xl’^{9)}o<i<n) the discrete-time filtration associated 
with X'^’- (resp. X^’-) 

:= a(X°;-, j<z) 

(resp. X/’-(6l) := a{Xl^^{e), j < i)) . 


3.2 Discrete-time scheme for (F, Z, U) 

We introduce an approximation of (F, Z) based on the discretization of (F°, Z^) and (F^, Z^). 
To this end we introduce the backward implicit schemes on tt associated with the BSDEs (12.41) 
and (1^ . Since the system is recursively coupled, we first introduce the scheme associated with 
(12.41) . We then use it to define the scheme associated with ()2.6I) . 

• Backward Euler scheme for (F^, Z^). We consider the implicit scheme (F^’-, F^’-) defined by 


[zir-.inm) 


giX^r^inm , 

[Y,Y{7:{e))AWf] , 7r(0) < , 


1 < z < rz , 


where Ej^’® = E[ . |Xji^’-(s)] for 0 < z < rz and s G [0,T]. 


(7r(d)),0)Af-, 

(3.3) 


• Backward Euler scheme for (F°,F°). Since the generator of (12.61) involves the process 
{Yt\t)) tG[ 0 T]! consider a discretization based on F^’-. We therefore consider the scheme 
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(yO.TT^ _2-o,7r) (defined by 


= g(X°n , 

^°-x = ^ !<*<«: 

where E° = E[ . \ for 0 < i < n, and P is defined by 

r{t,x,y,z) = f{t,x,y,z,Y^^^^{TT{t))-y), 

for all {t, X, y, z) € [0, T] X R X K X K. 

We then consider the following scheme for the solution (F, Z, U) of the BSDE (12.21) 

F*" = YX,tt<r + YX{7:ir))l,>r, 




(3.4) 


(3.5) 


for t G [0,T]. 


4 Convergence of the forward scheme 

We introduce the following assumption, which will be used to control the error between X and 

(HFD) There exists a constant K such that the functions b, a and /3 satisfy 

|6(t, x) — 6(t',a;)| + |(T(t,a;) — cr(t',a;)| < — 

|,5(t, x) —/3(<', x) I + |(T(t, x) — (T(t', x) I < K\t — t'\, 

for all {t,t',x) G [0,T] x [0,T] x R. 

In the following we provide an error estimate of the approximation schemes for X^ and X^ 
which are used to control the error between X and X'^. 


4.1 Error estimates for and 


Under (HF) and (HFD), the upper bound of the error between X'^ and its Euler scheme X^’^ 
is well understood, see e.g. [16], and we have 


E 


I V^O.TT I 2 

sup 


< K\tt\ , 


(4.1) 


foGlO.T] 

for some constant K which does not depend on tt. 

The next result provides an upper bound for the error between X^ and its Euler scheme X^’’" 
defined by (13.11) . 

Theorem 4.1. Under (HF) and (HFD), we have the following estimate 
sup e[ sup \xl{e) - Xl'^X‘n{e))\^] < K\tt\, 

eG[o,T] foG[e.T] 


for a constant K which does not depend on tt. 
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Proof. Fix 0 S [0,T], we then have 


E 


sup \xl{0) - 
'-t^[e,T] 


< 2E 
+ 2 E 


sup \xl{e) -xl{i,{e))[ 

'-t&[e,T] 


sup \xl{n{e))-xl';^{TT{e))\^ . (4.2) 


We study separately the two terms of the right hand side. 

Since 7r(0) < 9 < t, we have by definition XI{tt{0)) = X^ for all s € [O,7r(0)), and Xl{9) = X^ 
for all s G [0, 9), which implies 


X^i9)-X^i7ri9)) = 


/ b{s,X^,)ds 

J 7t(6) 


T(e) 

+ f cj{s 
Je 

-f 

Jtv{ 9 ) 




{s,X°)dW,+p{9,X°)+ / b{s,Xl{9))ds 


-/3(^(0),XO(,)) - f b{s,Xl{7ri9)))ds 

J-K(e) 

a{s,Xl{^{9)))dWs , 


for all t G [0, T], 

Hence, there exists a constant K such that 


\X^{9)-X^{n{9)f < k\ 



2 

/ b{s,X°)ds 
Jtv{B) 


‘ b{s,Xl{9))-b{ 

+ 

l-B 

/ a{s,X°)dWs 


l^iB) 


(a{s,Xl{9))- 


b{s,Xl{Tr{9)))ds 


ds 

a{s,Xl{TT{9)))dWs 
dWs^ 


T(e) 

'MO] 

f 

UiB) 

+ |/3(0,XO)-/3(^(0),XO(,))|"} . 

From (HF) and (HFD), we have 

E\M0,X°B)-PiA0)M°(9))f < if(kP+E|xO-X°(,)|") . 

We have from (HF) and (1^ 


(4.3) 


E 


r6 2 

^0' 


iTTiB) 


b{s, Xg)ds 


+ 


r® 2n 

"0\ 


Itv{B) 


a(s,X,")dW, < K\Tr\ 


which implies in particular E|Xg — < K\'k\ and hence 

E|/?(0,X°)-^(7r(0),X°(,))|2<iF|7r|. 
We have also from (HF) and (I2.10|) 


E 


iTT{e) 


h{s,Xl{-K{9)))ds 


+ 


l-niB) 


r(s,X](7r(6»)))dWs < K\-k\ 


Combining these inequalities with (14.3p . (HF) and BDG-inequality, we get 


E 


sup \Xli0) - Xlin{0))[ 

ue[B,i\ 


< K 


( fE\ sup |Xi(0)-Xi(^(0))|"lds + k|) . 
^Jb UG[e.s] J ^ 
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Applying Gronwall’s lemma, we get 


E 


sup |Ai(0)-Ai(7r(0))|" < i^|7r|. 


(4.4) 


To find an upper bound for the term E[sup(g[g|A(^(7r(0)) — A^(j^(7r(0))p] we introduce the 
scheme A’^(7r(0)) defined by 

r KieMO)) = KieMO)) , 

\ Xl{iT{e)) = X-_^(^(0))+6(^._l,A-_^(7^(0)))A^- + a(<._l,A-_^(7^(0)))AW,^ U>it{0). 
We have the inequality 


E 


sup |Ai(7r(6»)) - A^(j)(7r(6»))|^ 
Geie.T] ^ ’ 


< 2E sup \xl{n{e))-x:^,^{i,{e))\ 

Lte[e,T] 

+ 2e[ sup \x:^,^{7t{ 9)) - xl’;^i7T{e))\^] . (4.5) 
Ge[e.T] J 


Since X’^(7r(0)) is the Euler scheme of A^(7r(0)) on [7r(0),T], we have under (HF) and (HFD) 
(see e.g. [TB] ! 


E 


sup |Ai(^(0))-X;(,)(7r(0))|"j < i^(l+Ef|Ai(,)(7r(0))|"j)|7 

te[e,T] J ^ L J/ 


for some constant K which neither depends on tt nor on 0. From (12.101) . we get 

I 2 ' 


E 


sup |Ai(7r(0))-X;(,)(^(0))|" < K\n\ 

'-tefe.T] J 


(4.6) 


for all e e [o,r]. 

We now study the term E[sup(g[g j.] |A()'^jj(7r(0)) — A^(j^(7r(0))|^]. We first notice that we 
have the following identity 


E 


sup |A;(t)(7r(6»)) - A^(j)(7r(6»))p 


= E 


sup |A;(t)(7r(6l)) - A^(j)(7r(6»))|' 

'-tG[7r(e).T] 


Hence we can work with the second term. From the definition of X"^ and X^’’^, we get 


sup |A;(„)(7r(6»)) - xI'^^^{tt{9))\^ < 

ttG[7r(e),t] ^ ' 


rTT{t) 

K{\Xl^,){7r{9)) - Xlll^{n{9))\^ + / 5(^(s), A;(,)(^(0))) - 6(^(s), Aij)(7r(0))) 

J ■7r(0) 


ds 


sup 

UG['7T{9),t] 


^7r(u) 


J-K{e) 


(a{7T{s),X^^s)i^{9))) - a{Tr{s),Xl’^l^{TT{9)))'jdWs ) . 


Then, using (HF) and BDG-inequality, we get 


E 


sup \x:^^){n{9)) - Xl’i:^^{7r{9))\^] < K(E\xl^,){n{9)) - Xll-^i7T{9))\‘ 

ug[7r(e),t] ' 


E 


7r(0) ^uG[7r(0),s] 


sup |A;(^)(7r(6»)) - A^(^)(7r(6l))|^jds) . 
\^(a\ cl ^ ^ ^ -I / 


From Lipschitz property of /3, we have 


= E|X“(,) +^(7r(0),A°(,)) - -^(7r(0),A°(;))|^ 
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This last inequality with dUI) gives 


< K\n\. 

Applying Gronwall’s lemma, we get 
E 


sup A;i(t)(7r(6»)) - A^’(^)(7r(6l)) < K\t:\ . 

'-te[7r(e).T] J 


Combining (I4.2L (14.41) . (14.5|) . (14.61) and (14.7L we get the result. 


(4.7) 

□ 


4.2 Error estimate for the FSDE with a jump 


We are now able to provide an estimate of the error approximation of the process X by its 
scheme X'^ defined by (13.21) . 

Theorem 4.2. Under (HF) and (HFD), we have the following estimate 

I 2 ' 


E 


sup \Xt-Xf 

'-tgfO.T] 


< K\'k\ 


for a constant K which does not depend on tt. 

Proof. From the definition of (DH) and m we have 


E 

sup |At-A,-|" 

< E 

sup ^ 7 r(t) 

+ E 

sup A1 (t) - A^’(^)(7r(T)) ^ 


CgIo.t] 


Lie[o,T) 


Cg[t.T] J 


< E 


I vO X^O.TT I 2 

sup I A, -A^ ,)| 
Celo.T] 


f E^ 

0 Ce[6»,T]' 


sup |Ai(0)-Aij)(7r(0))|V(0) 


d9 


< A(|7r|+ sup e[ sup |a]( 6I) - A^(^)(7r(6l))|^j) . 

^ 6»e[0,T] ‘-se[6».T] 


From Theorem 14.11 we get 


E 


sup |Ai-A,- 

tg[0.T] 


< A|7r| . 


□ 

5 Convergence of the backward scheme in the Lipschitz 
case 

To provide error estimates for the Euler scheme of the BSDE, we need an additional regularity 
property for the coefficients g and /. We then introduce the following assumption. 

(HELD) There exists a constant K such that the functions g and / satisfy 

\g{x) -g{x') \ + \f{t,x,y,z,u) - f{t\x',y,z,u)\ < A"(|a; -a;'| + \t-t'\^) , 

for all {t, t', X, x', y, z, u) € [0, T]^ x x R x K x K. 

We are now ready to provide error estimates of the approximation schemes for and 

(V^,Z^), and then for (V,Z). 
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5.1 Regularity results 

In this part, we give some results on the regularity of the processes and We denote 
:= cr{X° , 0 < s < t} and := a{Xl{6) , 9 < s <t}. 


Proposition 5.1. Under (HF), (HFD), (HBL) and (HELD), there exists a constant K 
such that 


E 





< K 


(l+E[|xi(0)|y)| 


(5.1) 


for all 9 G TT. 


Proof. We first suppose that b, a, f and g are in C^. Let us define the processes A and M by 


and 


At := exp(^y dyf{Ql{9))dr'^ , 


Mt := 1+ [ Mrd,f{Ql{9))dWr , 

Jo 


where 0^(9) := (r, Xf(9),Yf^(9), Z^(9), 0). We give classically the link between V^Xf(9)(:= 
dXf(9)/dXg(9)) and (DsXf(9))ff<s<t the Malliavin derivative of Xf(9). Recall that X^(9) 
satisfies 

Xl(9) = X^{9)+ fb{r,X^(9))dr+ fa{r,X^{9))dWr , 9<t<T. 

Je Je 


Therefore, we get 

V^Xl{9) = 1+ [ d^b{r,Xl{9))X<>Xl{9)dr 
Jo 


d,a{r,X^i9))y^X^i9)dWr , 9<t<T, 


and for 0 < s < t 


DsXl{9) = a{s,Xli9))+ f d,b{r, XX9))D,X^i9)dr + f d,a{r, XX9))D,X^i9)dWr . 

J S J S 


Thus, we have 


DsXl{9) = V^Xi(0)[V®Xl(0)] V(s,Xi(0)). 


(5.2) 


Using Malliavin calculus we obtain that a version of Z^{9) is given by {DtYt^{9))^^^g j,y By Ito’s 
formula, we get 


AtMtZl{9) = E 


MT(ArVg(Xi(0))A^T(^)+y a./(0j(0)) A^'(0)A,dr) | J-i(0)] , 

for t e [9,T]. Using (15.2p . we get 

AtMtZl{9) = E[MT(ATVg(Xi(0))V®Xi(0) + y U,A,dr) | J-i(0)] [V®Xi(0)] ■^(t, Xi(0)) , 
with Fr ■= da:f{&l{9)^X^Xf{9). This implies that 

AtXItZl{9) = {^[G\Fl{9)]-Mt j FrArdr')[A7^Xl{9)\-\{t,Xl{9)) , 

J 0 
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with G := g{X^{9))y^X\,{0) + FrArdr^. Since b, a, f and g have bounded 

derivatives, we have 

E[|Gr] < oo , p > 2 . (5-3) 

Define := E[G|J'^(0)] for r € [9-,T]. From (15.31) and Doob’s inequality, we have 


ll"i||sp[e.T] < oo , p>2 . 

Hence, there exists a process (j) such that 

m, = E[G\F^{e)] + [ (jjudWu , r € [6, T], 

Je 


and 


(5.4) 


^xUe)] ' 


||</>||ffp[e,T] < oo , p>2. 

We define Z by 

Zt{e) := {rnt - Mt [V® 

By Ito’s formula, we can write 

Zt{0) = Ze{e) + [ alie)dr + [ al{0)dWr , 0 <r <T . 

Je Je 

Since 6, cr, / and g have bounded derivatives, we get from (|5.4p 

sup \\Z{0)\\^g^,gj., < CX) , P>2, 

6»e[0,T] 

and 

sup ^ ^||a^( 0 )||//p[e_ 7 ’] + ||Q:^(6*)||_f/p[6i,T]) < oo , p >2 . 

We now write for t € [ti, ti+i) 

E[\Zt{0) - Z,M\"] < + 

with 

r II, := E[\Zt{0) - ZtM\"\a{U,Xl{0))\^] , 

\ll, := E[\Ztm"\a{t,Xl{0)) - a{U,Xl{0))\^] . 

We give an upper bound for each term. 


lit = E 


< KE 


E[\Zti0) - ZtM\"\Tli0)]\a{U,Xli0))f 

r 

{\al{0)\'^ + \al{0)\'^)dr sup \a{t,Xl{0))\^ 
i te[0,T] 


which implies 


pii+i p r 

Il^^dt < K\-k\E / {\al.{0)f + \almldr sup \a{t,Xli0))\^ 

ti ^Jti te[9,T] 


(5.5) 


(5.6) 
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therefore we have 


< K\Tr\E / {\al{e)\^+ \al{0)\^)dr sup |cr(t, | 


'-Je 


1=0, ti>0 '' 

From Holder’s inequality and (HFD) and (HF), we have 


te[6»,T] 


E ^Ltdt < K\7 t\e\ r {\aU9)\'' + \c^l\\0))drY(l + E 

n * a ti J 6 


2=0, ti>e 

Using (EH), we get 


sup \xl{e)\ 

'-K[e,T] 


n—1 

E 

2=0, ti>e 


iLtdt < 


< 


K\tt\(i 


(l+E 

sup 1 

(l + E 

'\xl{0) 




We get from (15.51) . (HFD) and (HF) 


lit < K IE 


(E[\Zt-zl‘^\xl{e)\^] +E[\xl{e)Zt-xl{0)Ztf]+\TT\^^ . 


Arguing as above, we obtain 

"-1 rti+i 


E 

i=o, ti>e ‘ 


E 


\Zt-Ztr\Xl{0)\ dt < K\Tr\E sup (l + |At^(6»)H 
-I ^9<t<T 


Moreover, from Ito’s formula, X^{0)Z is a semimartingale of the form 

Xl{0)Zt = X^{0)Ze+ [ aldr+ f d^dWr , 

Je Je 

where ||o;^||//2[0 (7-] + ||o;^||j/2[0 < A"(l + E[|X^(^)|^]^). Therefore, we have 

E\\Xl{0)Zt-Xl{0)Zll < Ke\ [+ \al\'^)dr 


L 


(5.7) 


which implies 


n —1 

E 

2=0, ti>e^ 


E 


Xl{0)Zt - Xl{0)Zl^ < K\tt\E (1 + \xl{0t) 


(5.8) 


Using EH and (15.81) we get the result. 

When 6, tr, /3, / and g are not in , we can also prove the result by regularization. We first 
suppose that / and g are in Cl- We consider a density q which is on R with a compact 
support, and we define an approximation of {b,a,l3) in by 


(&^cr^/3'')(^,a;) = ^ Jj^b,a,l3){t,x')q( J^ ^dx , (t,a;)€[ 0 ,r] 


We then use the convergence of {X^''^{0),Y^’^{0), Z^'^{0)) to {X^{0),Y^{9), Z^{0)) and we get 
the result. Next we assume that / and g are not Cl and we consider for that and g^ which 
are defined as previously and we get the result. □ 

Using the link between and Ag(0), we obtain that the bound (15.11) is actually uniform 
in 0. 
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Corollary 5.1. Under (HF), (HFD), (HBL) and (HELD), there exists a constant K such 
that 


E 


zlis)-zl^,^{e)\^dt < k\7t\ 


(5.9) 


for all 9 G TT. 

Proof. Since is a Brownian diffusion, we have for any p > 2, from (HFD) and (HF), that 

E 


sup |X0|P 

'-tG[0,T] 


< OO . 


We notice that from the Lipschitz property of /3 we have 


E 


\XI{9)\ 


= E 




< k(i + e 


sup |X°| 
'-te[o.T] 


< OO . 


Combining this result with (EU, we get (EH) 
We now study the regularity of 


□ 


Proposition 5.2. Under (HF), (HFD), (HBL) and (HELD), there exists a constant K 
such that we have 


E 




I yO 7 O I 


'dt 


< K\tt\ . 


Proof. The proof is similar to the previous one. The only difference is that the BSDE (12.6|) 
involves Y^. We denote 0° = (r, Yf), Y)^{r) — !)?)■ We first suppose that b, a, (3, f and 
g are in . We recall that 


= g{X^^) + f{Ql)ds- 


Therefore, for 0 < r < t < T, we have 


DrY(> = Vg{X°)DrX°+j^^(d,f{e°)DrX° + {dy-du)f{e°)DrY° 

+dJ{e°)D,Z° + 5„/(0°)D,r/(s))dr - ^ DrZ^dWs , 

where DX(!, DY(), DZ^ and DY)^{r) denote the Malliavin derivatives of X^, 1)?, Z^ and Y)^{r) 
for r G [0, T]. Using Malliavin calculus, we obtain that a version of Z° is given by (Dty)°)tg[o.T]- 
By Ito’s formula, we get 


AtMtZt — E 


MT(^ATVg{X°)DtX°+ f (a,/(0°)AX° + 5„/(0°)Ab;i(r))A,dr 




dx" 


where A* := exp(f*(dy—du)f(0^)dr) and Alt ■= 1+/q Afrdzf(0^)dWr. Denote by := 
and VX((9) := for 0 < t < 9 < T. We then have for r < s < T 

DrX)(s) = (l + d,/3(s,X°))DrX° = (l + d,/3(s,X°))VX^a(r,X°)[VX°]-^ , 


thus we can see that DrX)(s) = VX^(s)CT(r,X°)[VAr°] Therefore, we get by writing the 
SDEs satisfied by (DrX)(9))s^iff^T] for r < 9, and (VAi^(0))sg[6/,T] 

DrX)(9) = VX)(9) [VX°] ■ V(r, X°) , r<9<s. 
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Writing the BSDEs satisfied by {Dj.Y^{9))s^]^g rp] for r < 0 and (VF/( 0 )) 5 g[g 7 -], and using the 
previous equality, we get 

DrY} (s) = vy/(s) [VXO] ■ V(r, , s < 0 . 


This implies 

^tMtZt = E 


MT(ATVg(X°) 


VXS+ / Fr-Ardr 


[VX°] V(t,X°), 


with Fr := 9a;/(0°)VX° + i9„/(0°)Vy^^(r). We can write 


AtMfZt — 


(e[G|J-°]-^ MtF.A.dr) [VX°] V(t,X°), 


with G := Mt{AtAI g{X^)VFrArdr). Since 6 , cr, / and g have bounded derivatives, 


we have 

E[|G|P] < 00, p>2. 

Define := E[G|X°] for r € [0,T]. From (15.101) and Doob’s inequality, we have 

ll"i||sj’[o,T] < 00 , p>2 . 

Hence, there exists a process </) such that 


and 


We define Z by 


— E[G] + f (j)udWu , r € [0,T] , 
Jo 


ll<^l|ff!>[o,T] < 00, p>2. 


Zt := {AtMt)-^(mt-Mt j F,A,dr) [VX^] ^ te[0,r]. 

By Ito’s formula, we can write 


Zt — Zq + f Q.^ds F f cY^dWr , t G [0, T'\ . 
Jo Jo 


Using the fact that 5, tr, / and g have bounded derivatives and (15.111) . we get 

lsP[ 0 .T] 

and 


II^IISprn.Ti < 00^ P> 2 , 


l|Q^^I|jrp[o,T] + l|o^llffp[o,T] < 00 , p>2. 

We now write for t £ [ti, ti+i) 

E[|ZO-y°j 2 ] < 

with 

r E[\Zt-Ztf\a{U,Xl)\^] , 

\lly.= E[\Zt\^\a{t,X?)-a{U,Xl)\^] . 


(5.10) 


(5.11) 


(5.12) 
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As previously we give an upper bound for each term. 


lit 


< KE 


sup 

te[o,T] 




From Holder’s inequality and Lipschitz property of tr, we have 


n—1 


E 



< 


Ar|7r|]E 




sup 

tG[0.T] 




Using (I5.12I1 . we get 



i+i 




< 


K\Tr\ . 


From (HFD) and (HF), we get 

+ kp) ■ 

Arguing as above, we obtain 


< K E 


\Zt-Zt. 




-E 


\X?Zt - X^ Z. 


n—1 

E 


i=0 



^t - Zt, 



< 



sup 

t6[0,T] 





Moreover, X^Z is a semimartingale of the form 


X°Zt = X°Zo+ [ aldr+ [ a^dWr 
Jo Jo 

where + ||Q:^||_f/2[Q yj < K and we have 


E 


Xl>Zt - Xl\Z. 


< KE 


^)dr 


which implies 


E 


i=0 

When 6, a, f and g are not C^, 

EH 


E 


XJZf-XJ.Zt. 


< K\tt\ . 


we can also prove the result by regularization as for Proposition 

□ 


5.2 Error estimates for the recursive system of BSDEs 

We first state an estimate of the approximation error for (Y^, Z^). 

Proposition 5.3. Under (HF), (HFD), (HBL) and (HELD), we have the following esti¬ 
mate 


sup 

eG[0,T] 



^7r(t) 


(k{9)) 


2 



for some constant K which does not depend on tt. 




< A:|7r| , 
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Proof. Fix 0 S [0,T] and t G [0,T]. We then have 


E 




< 2E[|y/(0)-y,i(7r(0))| 

+ 2E[|y/K0))-r\’f)K0))| 


(5.13) 


We study separately the two terms of right hand side. 

Define 5Xl{0) := Xl{e) - Xl{n{e)), 5Y,\e) := r/(0) - r,i(7r(0)) and 5Zl{0) := Zl{e) - 
ZI{'k{9)). Applying Ito’s formula, we get 




ds 


+ 2^ 5Y}{0)SZl{e)dWs+ \SZl{0)\^ds, 
where Ql{0) := {s,Xl{0),Y}{0), Zl{0),O). From (HBL) and (HELD), we get 

E[|<5F/(0)P] < k(e[\SXU0)\^] +e[I^^ \SY,H0)\\6Xl{0)\ds\ + e[ \SY,\0)\^ds 


+ E 


(0) I (0) Ids]) - E [ ^ \6Zl (d) pds 


Using the inequality 2ab < a^/rj + ■qb'^ for a,b gR and q > 0, we can see that 

E[\SY,\0)\^] +e[I^^ \SZl{0)\^ds] < K(E[\SXm\"] + ^^E[|dr/(d)|2]ds 


- E 


|dA;(d)|2ds]) . 


From (14.41) and Gronwall’s lemma, we get 

E[|r,i(d)-y/(^(d))|"] < K\7r\. 


(5.14) 


(5.15) 


We now study the second term of the right hand side of (15.131) . Using the same argument 
as in the proof of Theorem 3.1 in [3], we get from the regularity of Z^ given by Corollary 15.II 

(5.16) 


E 


|y/w^))->;\’r)W^))r] < Ku\. 


This last inequality with (j5.13() and (I5.15|) gives 


sup < sup E 
6»e[o.T] 1 te[e,T] 


|r,i(d)-r4y7r(d))|"]} < KU\. 

We now turn to the error on the term Z^ (0). We first use the inequality 


E 


'-Je 


\Zli0) - 


< 2 E 


'-Je 


\Zl{n{0)) - Zl’;^{7r{0))\^dt 


■ 2 E 


\dZl{0)\^dt 


(5.17) 


Using (15.1411 and (15.151) with t = 0, we get 


E 


\5Zl{0)?ds 


< AIttI . 


(5.18) 
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The other term in the right hand side of (I5.17P is the classical error in an approximation of 
BSDE. Therefore, using Corollary IS.ll and (I5.16|) . we have 


E 


^Je 




< Kh\ . 


(5.19) 


Combining (15.171) . (15.181) and (15.191) . we get 


E 


Je 


zl{e)-zl^;^{n{e))\-‘dt < K\n\. 


□ 


We now turn to the estimation of the error between and its Euler scheme (13.4|) . 

Since this scheme involves the approximation Y^''^ of we first need to introduce an interme¬ 
diary scheme involving the ’’true” value of the process . We therefore consider the scheme 
(yo.TT^^o.TT) defined by 

Yt\ = + / {U-I. , Y,l\ , , Y,l^ ) At- , ^^ 20 ) 

[Yt^W-] , 1 < * < n . 


Using the regularity result of ProDOsition l5.2l and the same arguments as in the proof of Theorem 

3.1 in i, we get under (HF), (HFD), (HBL) and (HELD) 


sup E 
te[o,T] 


I vO _ 1^ 


-E 


■-Jo 


I 7 O _ /O.ir |2 j. 


< K\7t\ . 


(5.21) 


With this inequality, we get the following estimate for the error between (y°, Z^) and the Euler 
scheme (IXil) . 

Proposition 5.4. Under (HF), (HFD), (HBL) and (HELD), we have the following esti¬ 
mate 


sup E 
ie[0.T] 


0 w0,7r|2 


y" - y 


7r(t) I 


-l-E 


- z^fj^dt 


^JO 


7T{t) I 


< K\Tr\ , 


for some constant K which does not depend on tt. 
Proof. We first remark that 



0 T^0,7r|2 


yU _ 


7i-(i) I 


yO _ |2 j. 


< 2 sup E 

te[o.T] 

~T 

< 2E' 


I vO _ P 

Ut 


'-v'O 


7 O _ rti 


-I- 2 sup E 
te[o.T] 

• T 

2E' 


|-yr0,7r -^0,7r I 2 

r7r(t) ■^7r(t)| 


/O 


r^0,7r r^Oj'TT I ^ 

^TTp) “ ^7r(£)| 


Using (|5.2ip, we only need to study suptg[o,T] 'H\Y°'iJ^-Y°'-^\^] and E[/J" fZ^’-^ - Z°’-^\^df\. To 
this end, we need to introduce continuous schemes for all 0 < f < n — 1. Since E[|y(°’^p] < 00 
and < 00 for all 1 < z < n, we deduce, from the martingale representation theorem, 

that there exist square integrable processes Z®’’"' and such that 


y*y = E[Y°f-^\Fu] + Z°’-dW,, 

^ [YuZ I iTdWs . 


YC’- = E 
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We then define 


yO,. ^ yO,. _ Y°:^, Z°:^, Y,]’^{t,) - Y°’^) + J* Z°'^dWs , 


yo,. ^ yo,. _ _ t,)f{u,x^:-,Y,^:-,z°:\Y,]{u) - y°’^) + I ffdWs, 


for t G [ti, ti+i). Let i € {0,..., n — 1} be fixed, and set SYt := Y^’"^ — Y^’^ , 6Zi := Z^’^ — Z^’’^ , 
S^Z,^.= and 5ft := /(t., X°;-, Z°;^ Y^:%U)-Yl^-)-f{U, Z°;^ Y^iU)- 

Y^f^) for t € [ti,ti+i). By Ito’s formula, we compute that 


/ ti+l r'i + l 

E\5Zfds - ^ / E[6YJfs]ds , k<t< U+i . 

Let a > 0 be a constant to be chosen later on. From the Lipschitz property of / and the 
inequality 2ab < aaf + 6^/a, we get 


rti+i rti+i p 

At < a E\5Y,\'^ds + - E \dYtf + \SZ,f + \Yt]{U) - Y,]’^{U)\‘‘ 

Jt a Jt 1 


ds . 


Using Proposition 15.31 we get 


At < a 
We can write 


K K K 

E|<5y«pds + — |7r| E|,5Ftj2 _ / Ej^Z.pds + — Itt 

a Jt ^ 




E\5YtY < E\5Yt,^,Y+ J E\dZj^ds < a 


/ ti+i 

El^nl' 


ds “b Bt 


(5.22) 


where 


B, := E\5Yt,^,\'^ + —\7r\E\5Z,\^ + -\7r\E\5Ytf + —\7r\^. 

a a a 

By Gronwall’s lemma, this shows that E|(5ytp < for ti <t < ti+i, which plugged in the 

second inequality of (15.221) provides 

E\5Yt\^ + j^'*\\5Zfds < B,(l + a|7r|e“l’^l) . (5.23) 

Interpreting (resp. Z^^^) as the projection of (resp. in H'^[ti,ti+i] on the set of 

constant processes, we have 


/ ti+i 

E\5Z,f 


ds < 


ii+1 


E\5Zfds . 


(5.24) 


Applying (j5.23l) for t = U and a = 2K, and using the previous inequality, we get 

/ ti+i 

E\5Zfds < k2{^)E\5Yu^,\^ + h{^)\-K\’^ , 0<i<n-l, 


"'here h ('-t'I — 5, / x _ l+2if|7r|e^^l’'l , . / x _ ^+K\ 

Where fci(^) - m ^2)^) - i_ M _^|^|2,2K|,| and 

Since for small |7r| we have fci(7r) > 0, we get 


,2K|,x| 


_J^_/sr|7r|2e2K|’f| ' 


ElJUtpp < k2{TT)E\5Yu^,\^+ ks{TT)\^\^ , 0<i<n-l 
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for \'k\ small enough. 

Iterating this inequality, we get 

n 

3=i 


Since A:2(7r) > 1 and 6Yt^ = 0, we get for small |7r| 

E|5YtJ^ < |7r|fc3(7r)fc2(7'‘)T^ < A"|7r| , 0<i<n, (5.25) 


which gives 

sup E[K’r)-y:(r)r] < K\A. 

tG[0,T] L ^ ^ ^ ^ J 

Summing up the inequality (15.231) with t = ti and a = 2K and using (15.241) . we get 

1 p T ^ 1 

< 2K\7r\e^^\^\J2E\6Ytf + {l + 2K\7r\)E\SYt„ 

I n—1 

+ (- + iTTle^^H) (kI + K| ^ E\SYu p) . 

Using (|5.25p . we get for |7r| small enough 




nzii:^-K[:/ds < k\a 


□ 


5.3 Error estimate for the BSDE with a jump 

We now give an error estimate of the approximation scheme for the BSDE with a jump. 

Theorem 5.1. Under (HF), (HFD), (HBL) and (HBLD), we have the following error 
estimate for the approximation scheme 


sup E 

\Yt-Yf^\^ 

+ E 

pT 

/ \Zt-Zf\^dt 

+ E 

pT 

/ \t\ut-u:?dt 

tG[0.T] 



^Jo ^ 




for some constant K which does not depend on tt. 


Proof. 

Step 1. Error for the variable Y. Fix t € [0,T]. From Theorem 12.II and (13.5p . we have 


E 


Yt-Y^ 


= E 


IwO _ I 
Ut 


'^t<T 


-E 


|F/(r)-y\’;;(7r(r))pi 


t>r 


Using (DH), we get 


E 


Y-Yf' 


< E 


Y}>-Y' 


0,7r I 2 

(t)l 


E 


r/(0)-r\’")(7r(0))|\*>,7T(0) 


< it: E 


Y^-Y' 


0,7r 12 


Tv{t) \ 


+ sup sup E 
eG[o,T] se[e,T] 


\Y,\O)-Y^i:^{n{0)) 


Using Propositions 15.31 and 15.41 and since t is arbitrary chosen in [0,T], we get 


sup 

te[o.T] 


E 


\Y-Y,^ 


< K\ti\ . 
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step 2. Error estimate for the variable Z. From Theorem 12.11 and (I3.5L we have 



T F 9 1 




r r'^ , 2 1 

E 

/ \Zt-Zf\Ut 
^Jo 

= E 

1 \z?-zl';^\dt 

+ E 

/ \zU'r) - dt 

^JTAr 


Using (DH), we get 


E 


\Zt - Zfl dt 


'-Jo 


ff. 

pT pT 


I t Tr(t)\ 


7t(0) 


dtdO 


E 


< K 


10 J0 




Zl{0) - 7t(0) 


dtdO . 


I ryO ryOiTT I 

\^t ^7T{t)\ 


'dt 


From Propositions 15.31 and 15.41 we get 

pT 


E 


\Zt - Zf\ dt 


'-Jo 


sup E 
0e[o,T] 


< KU\ . 


^{0) - Z^(j)(7r(6»))|^ dt^ . 


Step 3. Error estimate for the variable U. From Theorem 12.II and (13.511 . we have 


E 


Ut-Un Xtdt 


'-Jo 

Using (HBI), we get 


< KE 


(|U/(t) - + \Yt° - YX/)Xtdt 


E 


Ut-ufrXfdt 


<K( sup sup E \Y^{0) — Yb'^J'jT{0))\ + sup E 

eefo.T] tG[6i,T] L ^ J te[o,T] 


Ut 


Combining this last inequality with Propositions 15.31 and 15.41 we get the result. 


□ 


Remark 5.1. Our decomposition approach allows us to suppose that the jump coefficient /3 is 
only Lipschitz continuous. We do not need to impose any regularity condition on f3 and any 
ellipticity assumption on /^ + V/3 as done in [4] in the case of Poissonian jumps independent of 

W. 


6 Convergence of the backward scheme for the quadratic 
case 

In this section we assume that (HBQ) holds and that a{t,x) = (j{t,0) = a{t) for any t G R+ 
and a; G R. 

Before giving the error of the scheme we give a uniform bound for the processes Z^ and Z^ 
which allows to prove that the BSDE (12.21) is Lipschitz and thus we can use Theorem 15.11 For 
that we introduce the BMO-martingales class, and we also give some bounds for the processes 
X°, X\ Y° and pi. 

6.1 BMO property for the solution of the BSDE 

To obtain a uniform bound for the processes Z^ and Z^ we need the following assumption. 
(HBQD) There exists a constant K / such that the function / satisfies 

|/(t, X, y, z, u) - /(<', x', y', z\ m')I < Kf [\x - x'\ + \y - y'\ + \u - u'\ + \t - t'\^] 

+LfAl + \z\ + \z'\)\z-z'\, 
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for all (f, X, x', y, y', z\ u, u') G [0, T] x x x x 


In the sequel of this section, the space of BMO martingales plays a key role for the a priori 
estimates of processes and . We refer to [T3] for the theory of BMO martingales. Here, 
we just give the definition of a BMO martingale and recall a property that we use in the sequel. 


Definition 6.1. A process M is said to be a BMOw[0,T]-martingale if M is a square integrable 
¥-martingale s.t. 


\m\ 

BMOw[0,T] 


sup E 

rGTifO.T] 


Mx — Mt 



< oo , 


where 7i’[0,r] denotes the set of¥-stopping times valued in [0,r]. 

The BMO condition provides a property on the Dolean-Dade exponential of the process M. 

Lemma 6.1. Let M be a BMOr[0,T\-martingale. Then the stochastic exponential £{M) de¬ 
fined by 

£{M)t = exp(Mi-i(M,M))^ , 0 <t<r, 

is a uniformly integrable ¥-martingale. 

We refer to m for the proof of this result. 


We first state a BMO property for the processes and which will be used in the sequel 
to provide an estimate for these processes. 

Lemma 6.2. Under (HF), (HBQ) and (HBQD), the Tnarthngales J'q Z^dW^ and J'q 
0 € [0,T] are BMOrlO,T]-martingales and there exists a constant K which is independent from 
9 such that 



/ Z^dWs 


< K, 



Jo 

BMOrlO,T] 



sup 

[ Zl{9)ls>edWs 


< K . 


6 le[ 0 .T] 

Jo 

BMOrlO,T] 



Proof. Define the function (f : 

IR ^ M by 




f{x) 


- oml , 

X G R . 

( 6 . 1 ) 


We notice that </> satisfies 

4>'{x) > 0 and -cj>"{x) — Kqcj)'{x) = 1 , 

for X > 0. Since and H^(.) are solutions to quadratic BSDEs with bounded terminal 
conditions, we get from Proposition 2.1 in m the existence of a constant m such that 

ll^°lls“[o.T] < m and sup ||ri( 6 »)|| 5 oo[e_T] < m. ( 6 . 2 ) 

eG[0,T] 

Applying Ito’s formula we get 


cl>{Y° + m) + e( £ + m)\Z°J^ds\F,) = 

E{f{Y° + + e( £ ^'{Y° + m)f{s, X°, Y°, Z^ - F°)ds| J-.) 
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for any F-stopping time v valued in [0,T]. From the growth assumption on the generator / in 
(HBD), (16.11) and (16.21) . we obtain 


+ m)+E(^J^ |Z°| 2 ds|jr,) < 

E((^(K0+m)|J-,)+E( r cj)'{Y° + m)Kg{l + 2\\Y°\\s^+ sup ||yi(0)|Uo=[e.T])ds • 

eG[0.T] ^ 


This last inequality and (16.21) imply that there exists a constant K which depends only on m, 
T and Kg such that for all F-stopping times z/ € [0, T] 


E 




< K . 


For the process Z^, we use the same technics. Let us fix 0 S [0,T]. Applying Ito’s fomula we 
get 


cl>{YXe{0)+m)+E(^£^^ + m)|Zi(0)pds| = 


O nT 

cj,\YX0)+m)fis,Xlie)Xi0),Zli0),O)ds J-.ve) 
!/ve ^ 


for any F-stopping time v valued in [0,r]. From the growth assumption on the generator / in 
(HBQ), dSU) and (EJ), we obtain 


HYX9{0)+m)+E( f iZlm^ds T,) < Ei^iY^{0)+m)\X,) 

+ e( r cj,'{Y}{e)+m)Kg{l + \\Y\0)\\s^[e^T])ds X,) . 


This last inequality and (16.21) imply that there exists a constant K which depends only on m, 
T and Kg, such that for all F-stopping times ly valued in [0,r] 


E 


\Zl{0)\Hs>eds 


T, 


< K . 


□ 


6.2 Some bounds about and 


In this part, we give some bounds about the processes A° and which are used to get a 
uniform bound for the processes Z^ and 

Proposition 6.1. Suppose that (HF) holds. Then, we have 


|VX 0 | := 

and for any 0 € [0, T] we have 

\V^Xl{0)\ := 


dXf 


dx 


< e 


dXl{0) 


dxi[0) 


LaT 


< e 


LaT 


\VXl{0)\ := 


dXf{0) 


dx 


< (1 + Lae 


LaT 


0<t<T, (6.3) 

, 0 <t<T , (6.4) 

)e^“'^ , 0 <t<T . (6.5) 
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Proof. We first suppose that b and P are w.r.t. x. By definition we have 

= 1+/ d^b{s,X°)\/X°ds , 0<t<T. 

Jo 

We get from Gronwall’s lemma 

|VX°| < , 0 < t < T . 

In the same way, we have 

X^Xl{0) = 1+ [ d^b{s,Xl{e))V^Xl{e)ds , 0<t<T, 

Je 


and from Gronwall’s lemma we get 




<t<T . 


Finally we prove the last inequality. By definition 

XXl{0) = 1+ [ d,b{s,Xl{0))XXl{0)ds + d,Pi0,X°)VX^ , 0<t<T. 

Jo 

Using the inequality (|6.3p . we get 

|VXi(0)| < 1 + Lae^^^ + [ La\XXl{0)\ds , 0<t<T, 

Jo 

from GronwalFs lemma we get 

\XXl{0)\ < (1 + , 0<t<T. 

When b and j3 are not differentiable, we can also prove the result by regularization. We consider 
a density q which is on R with a compact support, and we define an approximation (5*^, /3*^) 
of (6, /3) in Cl by 

{b\P^){t,x) = ^ J {b,P){t,x')q( J ^ 'jdx' , (t, x) e [0, T] x R . 

We then use the convergence of [X^'^,X^''^{0)) to {X°,X^{0)) and we get the result. □ 

6.3 Some bounds about and 

In this part, we give some bounds about the processes and which are used to get a 
uniform bound for the processes and Z^. 


Lemma 6.3. Suppose that (HF), (HBQ) and (HBQD) hold. Then, for any 0 G [0,T] 

|V®r/(6»)| := < e^^'^+^f'^^iKg+TKf) , 0<t<T. (6.6) 


dXl{0) 

Proof. We first suppose that 6, / and 5 are Cj w.r.t. x, yandz. In this case (X^(0), F^(0), Z^(0)) 
is also differentiable w.r.t. Xg{0) and we have 


VX\0) = VgiXl.{0))V<^X^{0)- l\'^Zl{0)dWs 


(6.7) 


V/(s, Xi (0 ), r/ (0 ), Zi (0) , 0) {X^Xl (0 ), V^Y} (0 ), V® Zl i0))ds, 
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for t G [9,T]. Define the process R{0) by 

Rt{e) := dyf {s,Xl{e),Y} (9), ZUe), 0)1 S>sds) , 0<t<T. 

Applying Ito’s formula, we get 

Rtie)V^Yt\e) = RT{e)Vg{X^{9))V<^X^{9) 

+ J^ Rs{e)d,f{s,xl{e)Xm,zl{9),0)v<^xl{9)ds 

- Rs{e)x^zl{e)dwl{e), e<t<T, 

where the process W"^{9) is defined by 

W}{e) -Wt - f dJis,Xlie),YXO),Zli0),O)ls>eds 
Jo 

for t G [0,T]. From (HBQD), there exists a constant K > 0 such that we have 

2 


( 6 . 8 ) 


(6.9) 


djis,xlie)Xid),zli9),o)t,>9dw, 


BMOrlO,T] 


(l+ sup E 

[ \Zl{9)\H,>gds 


^ i9e7i[0,T] 


-* 


< 


< 


K 


(l+ 1^ Zl{0)l,>edWs 


BMOwlO,T] 


< oo 


where the last inequality comes from Lemma 16.21 

Hence by Lemma l6H] the process f(/p dzf{s, A^(0), Y^{9), Z]{0), 0)1 s>gdWs) is a uniformly 
integrable martingale. Therefore, under the probability measure Q^(0) defined by 


dQ\0) 


dF 


:= £( ['dzfis,Xli0)Xid),Zli0),O)ls>edW,) , 0<t<T, 

Jo ■' t 


we can apply Girsanov’s theorem and (0) is a Brownian motion under the probability measure 
Q^(0). We then get from (16.81) 


Rt{0)X^Y; = EQi(e)[i?T(0)V5(Xi)V®Ai + RMd.fis,XlY^\zlO)X’^Xlds 

This last equality, (HBQD) and (16.41) give 

|V®T/(6()| < + TKf) , 0 <t<T . 


Tt 


( 6 . 10 ) 


When b, f and g are not , we can also prove the result by regularization as for Proposition 


□ 

Lemma 6.4. Suppose that (HF), (HBQ) and (HBQD) hold. Then, 

\^Y)^{t)\ < (1+ TTAT/) , 0<t<r. (6.11) 

Proof. Firstly, we suppose that b, /3, g and / are €( w.r.t. x, y and z. Then, for any t G [0,T] 

VL)i(t) := 




dx 


= xgix^{t))yx^{t) 


( 6 . 12 ) 


+ xf{s,xl{t), Y}{t),zl{t), o){vxl{t), vy/(t), XZl{t))ds 


XZl(t)dWs 
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Applying Ito’s formula, we get 

Rt{t)VYt\t) = RT{t)VgiX^it))VX^{t) 

+ ^ Rs{t)d,fis, Xl{t), Y}it),Zl{t), 0)VXlit)ds 

- Rs {t)yzl {t)dW^ (t) , 0<t<T , 

where the process W^{.) is defined in (16.91) . We have proved previously that W^{t) is a Brownian 
motion under the probability measure Q^(t). We then get 

Rt{t)VY,\t) = EQip)[i?r(<)Vg(Ai(<))VXi(t)+^ Rs{t)d,f{s,Xl(t)X{t),Zl{t),Q)VXl{t)ds 

This ISiSt iri6QU&lity, Riici (lOl) give 

|vy/(t)| < {l + Lae^-^)e^^-+^f^^{Kg+TKf) , 0<t<T. 

When b, f and g are not , we can also prove the result by regularization as for Proposition 

HH □ 

Lemma 6.5. Suppose that (HF), (HBQ) and (HBQD) hold. Then, 

|Vy 4 °| < e(2^^+^“)^(Xg+A:/r)(l+TA:/e^^^(l + Lae^“'^)) , o<t<r. 

Proof. We first suppose that b, (3, g and / are Cl w.r.t. x, y, z and u, then is 

differentiable w.r.t. x and we have 

VY° = Vg(A°)VX° 

+ r (v/(,.x.», z.«,y/(.) - y»)(vx.», vry vy.«. vy/,.) - vy").. 


Xt 


yz'idWs , 0 < t < T. 


Define the process by 


Rl := 


exp ( J\dy - a„)/(s, XO, r,o, ZO, y/(s) - Y°)ds) , o < t < t . 


Applying Ito’s fomula we have 

?0y7w 0 _ nO 


RlXYl’ = i?^Vg(X^)VX^ 

T 


+ ^ x°, ZO, (s) - Y°)XX^ds 

+ £ x°, zo, Y^is) - y,°)vr;(s)ds 


i?°VZ°dW° 


where dW^ := dWg —dzf{s, X°, y/(s) —Fg°)cis. From (HBQD), there exists a constant 

K > 0 such that we have 

2 

, ^ ^ , (q) 

s "> s ^ s ^ s 


dzfis, x°, r,°, Z°, Y^\s) - Y,°)dWs 


IQ 


< 


BMOr[0,T] 


h + sup E 

' r\z°\M 

x^> 

^ i3G7i[0.T] 


-* 


-Ml 


Z°dWs 


BMOr[0,T] 


< oo 


28 
















where the last inequality comes from Lemma 16.21 

Hence by Lemma EH] the process £(f^ dzf(s, Z°,Vj^(s) — Y^)dWs) is a uniformly 

integrable martingale. Therefore, under the probability measure defined by 



dF 


Tt 



a^/(s, x'i, Y^,z^, y/(s) - y^^ 



we can apply Girsanov’s theorem and is a Brownian motion under the probability measure 
Q°. Then, we get 




— Er 


R°TXg{X^)XX^ 


+ £ R°d,f{s, Xl yo, yo, F/(s) - Y°)VX°,ds 
+ R°d^f{s, X°, Y°, Z°, Y}{s) - Y°)VY,\s)ds 


Tf. 


Using inequalities (16.31) and (16.111) we get 


|Vrt°| < e^^^f+^‘^^'^{Kg + KfT){l + KfTe^f'^{l + Lae^<‘'^)) , 0<t<T. 


When 6, /3, / and g are not C^, we can also prove the result by regularization as for Proposition 

EH □ 


6.4 A uniform bound for and 

Using the previous bounds, we obtain a uniform bound for the processes and Z^. 

Proposition 6.2. Suppose that (HF), (HBQ) and (HBQD) hold. Then, for any 9 € [0,T], 
there exists a version of Z^ (9) such that 

\Zl{9)\ < + TKf)Ka , 9<t<T . 

Proof. Using Malliavin calculus, we have the classical representation of the process Z^{9) given 
by V^Y^{9)(y^X^{9))~^a{.) (see SectionEj)- In the case where b, f and g are Cl w.r.t. x, y 
and z, we obtain from (16.101) 

\Zlm < e^^^‘‘+^f'>^{Kg+TKf)Ka a.s. 

since |(V®X^(6)))“^ | < (the proof of this inequality is similar to the one of (16.41) 1. 

When b, f and g are not differentiable, we can also prove the result by a standard approximation 
and stability results for BSDEs with linear growth. □ 


Proposition 6.3. Suppose that (HF), (HBQ) and (HBQD) hold. Then, there exists a 
version of Z^ such that 

\ZQ\ < e‘^^^f+^-^^{Kg + KfT){l+TKfe^f^{l + Lae^‘^^))Ka, 0<t<T. 

Proof. Thanks to the Malliavin calculus, it is classical to show that a version of Z^ is given by 
Vy°(VX°)“^o-(.) (see SectionEj). So, in the case where b, /3, g and / are €( w.r.t. x, y, z and 
u, we obtain from (16.31) and Lemma 16.51 

\Z°\ < e2(^^+^“)'^(A'g + iLjT)(l + ri^/e^^^(l + Lae^“'^))M<, a.s. 

since |(VX°)“^| < (the proof of this inequality is similar to the one of (16.31) 1. 

When b, /3, g and / are not differentiable, we can also prove the result by a standard 
approximation and stability results for BSDEs with linear growth. □ 
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6.5 Convergence of the scheme for the BSDE 

Theorem 6.1. Under (HF), (HFD), (HBQ) and (HBQD) we have the following estimate 


sup E 

]Yt-Yf^f 

-l-E 

[ \Zt-Zf\^dt 

-l-E 

pT 

/ Xt\Ut-Uf\'dt 

te[o,T] 



^Jo ^ 


L Jo 


for a constant K which does not depend on tt. 


Proof. Fix M G R such that 

M > + Tii:/)iFa ; 

+ KfT){l + TKfC^f^il + , 

and define the function / by 


f{t,x,y,z,u) 


f{t,x,y,(pM{z),u) , 


{t,x,y,z,u) € [0, T] X R X M X R X R , 


where 


iPAliz) 


z if |z| < M 
if |2| > M 


We notice that cpM is Lipschitz continuous and bounded. Therefore we obtain from (HBQD) 
that / is Lipschitz continuous. 

Moreover, using Propositions 16.21 and 16.31 we get that under (HF), (HBQ) and (HBQD), 
{X, Y, Z) is also solution to the Lipschitz FBSDE 


Xt 

Yt 


x+ f b{s,Xs)ds+ f (T{s,Xs)dWs+ f j5{s,Xs-)dHs , 0<t<T 

Jo Jo Jo 

giXr) + ^ fis,Xs, Ys, f/,(l - H,))ds 


^ ZsdWs - [ UsdHs , 0 < t < T . 


Applying Theorem 15.11 we get the result. 


□ 
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